#include <bits/stdc++.h>
#define endl '\n'
using namespace std;
const long long mod=998244353;
long long n,m,q,k,ans=0,a[3000005],fact[3000005],inv[3000005],f[3000005];
long long pow2(long long a,long long b)
{
	long long ans=1,base=a;
	while(b>0)
	{
		if(b&1)
		{
			ans=ans*base%mod; 
		}
		base=base*base%mod;
		b>>=1;
	}
	return ans;
}
void init()
{
	fact[0]=1;
	for(long long i=1;i<=3000000;i++)
	{
		fact[i]=fact[i-1]*i%mod;
	}
	inv[3000000]=pow2(fact[3000000],mod-2);
	for(long long i=2999999;i>=0;i--)
	{
		inv[i]=inv[i+1]*(i+1)%mod;
	}
}
long long C(long long n,long long m)
{
	if(n<m)
	{
		return 0;
	}
	return fact[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	init();
	cin>>n>>m>>q>>k;
	if(n<m)
	{
		cout<<0<<endl;
		return 0;
	}
	for(long long i=1;i<=m;i++)
	{
		cin>>a[i];
	} 
	for(long long i=0;i<=m;i++)
	{
		f[i]=pow2(n-i,q);
	}
	for(long long i=0;i<=m;i++)
	{
		ans=(ans+(i&1?-1:1)*C(m,i)*f[i]%mod+mod)%mod;
	}
	cout<<ans*(n-m+1)%mod*(m<=q?pow2(k,q-m):0)%mod<<endl;
	return 0;
}
